Theory for large amplitude electrostatic ion shocks in quantum plasmas 
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We present a generalized nonlinear theory for large amplitude electrostatic (ES) ion shocks in 
collisional quantum plasmas composed of mildly coupled degenerate electron fluids of arbitrary 
degeneracy and non-degenerate strongly correlated ion fluids with arbitrary atomic number. For our 
purposes, we use the inertialess electron momentum equation including the electrostatic, pressure 
gradient and relevant quantum forces, as well as a generalized viscoelastic momentum (GVEM) 
equation for strongly correlated non-degenerate ions. The ion continuity equation, in the quasi- 
neutral approximation, then closes our nonlinear system of equations. When the electric field is 
eliminated from the GVEM equation by using the inertialess electron momentum equation, we then 
obtain a generalized GVEM and the ion continuity equations exhibiting nonlinear couplings between 
the ion number density and the ion fluid velocity. The pair of nonlinear equations is numerically 
solved to study the dynamics of arbitrary large amplitudes planar and non-planar ES shocks arising 
from the balance between harmonic generation nonlinearities and the ion fluid viscosity for a wide 
range of the plasma mass-density and the ion atomic-number that are relevant for the cores of giant 
planets (viz. the Jupiter) and compact stars (viz. white dwarfs). Our numerical results reveal that 
the ES shock density profiles strongly depend on the plasma number density and composition (the 
atomic-number) parameters. Furthermore, density perturbations propagate with Mach numbers 
which significantly dependent on the studied plasma fractional parameters, ft is concluded that the 
dynamics of the ES shocks in the super-dense degenerate plasma is quite different in the core of a 
white dwarf star from that in the lower density crust region. 



PACS numbers: 52.27.Gr,52.30.Ex, 52.35.-g, 52.35.Fp, 52.35.Tc 
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INTRODUCTION 



Recently, there has been a great deal of interest [lH8| in studying various collective phenomena in dense quantum 
plasmas which are ubiquitous in a wide-range of phy sical systems, including compact astrophysical objects (e.g. the 
cores of white dwarfs 0, [l(| and neutron stars jHf), giant planetary cores ( e.g. the interior of the planet Jupiter 
[Hj], high-energy density com pres sed plasmas created by powerful laser beams 13 1, as well as semiconductors (l4Lll5| 
and metallic nanostructures [16( for high-technological applications, table-top quantum free-electron lasers (FEL) 
fl7j ] for producing coherent tunable x-rays and gamma rays, and quantum diodes [ljj . etc. Quantum plasmas arc 
usually composed of degenerate electron and nondegenerate ion fluids. Since, in quantum plasmas the number 
density of degenerated electrons, no, is rather high, the average inter-electron/ion distances d ~ (3/47rno) 1//3 can be 
comparable with the electron thermal de Broglie wavelength, 2nh/ \/ m e Vr , which is much smaller than the Landau 
length L — e 2 /ksT e . Here where e is the magnitude of the electron charge, ks the Boltzmann constant, T e the thermal 
electron temperature, H the Planck constant divided by 2n, m e the rest mass of the electron, and Vt = {huTe/me) 1 ! 2 
the electron thermal speed. In quantum plasmas, the wavefunctions of degenerated fermionic ingredients (electrons) 
overlap giving rise to the Fermi-Dirac distribution function and density of states (DoS). In fact, the importance 
of the quantum electron degeneracy was recognized eighty years ago by Fowler [l9| and Chandrasekhar Q in the 
study of non-relativistic and ultra-relativistic pressures for degenerate electrons that were required for setting up 
the hydrostatic stability of a gravitating star that eventually becomes a compact star called a white dwarf. Clearly, 
in quantum plasmas, the equations of state (EoS) significantly differ from classical plasmas where non-degenerate 
electrons obey the ordinary Maxwell-Boltzmann distribution function. In quantum plasmas, the electron and ion 
gamma factors (which are defined as the ratio of the Coulomb potential to the thermal energy of the plasma particles) 
are given as T e = e 2 /dKBT e = L/d and T; = Z 2 e 2 / dksTi, respectively, being smaller and much larger than unity, 
implying a weak correlation for degenerate electrons and strong correlations for non-degenerate ion fluids in the 
quasineutral plasma regime. 

The collective behavior of a quantum Fermi plasmas with a statistical ensemble of degenerate electrons has been 
investigated in the pioneering papers of Klimontovich and Silin (20| , and later by Bohm and Pines [2l| , who have used 
the quantum kinetic Wigner and Poisson equations to obtain the dispersion relation for quantum electron-plasma 
oscillations (EPOs). The frequency spectra of the latter can be shown to be, w = (uj 2 e + k 2 Vp e + Ti 2 k A / 'Arrie) 1 / 2 , where 
u} pe — (47rrioe 2 /m e ) 1 / 2 is the electron plasma frequency and Vf b = (fi/m e )(37r 2 no) 1/ ' 3 is the Fermi electron speed. 
The frequency spectra of quantum EPOs have also been derived [l|, [22[ from the quantum electron fluid equations 
[23| , composed of the electron continuity, non-relativistic electron momentum equation including the quantum statis- 
tical pressure and the quantum recoil effect, and Poisson's equation. On the other hand, in his classic experimental 
work, Glenzer et al. [1 3j j have conclusively demonstrated the quantum mechanical collective effect by measuring the 
scattering of electromagnetic waves off the quantum EPOs in a high-energy density (warm dense matter) laboratory 
plasma. Numerical analyses of Poisson's equation and the nonlinear Schrodinger equation for non-relativistic degen- 
erate electron fluids revealed the formation of a dark soliton and two-dimensional quantized vortices Q, as well as 
nano-structures [24j in quantum electron plasmas with immobile ions. 

On the other hand, inclusion of the ion dynamics leads to electrostatic oscillations (ESOs) supported by the restoring 
force due to the quantum statistical pressure, the quantum Bohm force and the inertia of non-degenerate ions. By 
using Poisson's equation, the momentum equation for non-relativistic degenerate inertialess electrons and the classical 
hydrodynamic equations (the continuity and momentum equations) for non-correlated ions, Hass et al 25[ derived 
a modified Korteweg-de Vries (KdV) equation for small amplitude nonlinear ion oscillations in an unmagnetized 
quantum plasma. The KdV equation admits a localized solitary pulse as a stationary solution. Large amplitude 



localized ion oscillations in a collisionless quantum plasma were investigated by many authors }26H28| ignoring the 
ion correlation effects. The latter and the ion fluid viscosity have recently been incorporated in the small amplitude 
shock wave theory [29( of ion oscillations in a quantum plasma. 

In this paper, we present a theory for fully nonlinear nonstationary ES ion shock waves in a quantum viscous plasma 
with mildly coupled degenerate electrons and strongly correlated ion fluids for a wide-range of the Chandrasekhar's 
relativity parameter and the plasma atomic number. For our purposes, we use the hydrodynamic force balance 
(the electric force balancing the gradient of the degenerate electron pressure, gradient of the potentials involving 
the electron-exchange and electron-correlations, the gradient of the Bohm potential) for inertialess electrons and the 
generalized viscoelastic ion momentum (GVIM) equation, together with the ion continuity equation, invoking the 
quasi-neutrality approximation. The GVIM equation contains the essential physics of the ion correlation decay rate, 
as well as of the ion fluid bulk and shear viscosities. When the electric force is eliminated from the GVIM equation 
by using the electric field, from the momentum equation for inertialess degenerate electron fluids, we obtain a pair of 
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equations exhibiting the coupling between the ion fluid velocity and the ion number density. The pair of equations 
are then numerically solved to study the dynamics of IS ion shocks in planar as well as non-planar geometries. Our 
numerical analysis exhibits distinguished features for the formation of ES shock structures in different ion mass-density 
ranges which might play a very important role in acceleration of charged particles in dense quantum plasmas, such 
as those in the cores of white dwarf stars and in compressed plasmas created by intense laser and electron beams. 



THE VISCOELASTIC HYDRODYNAMICS MODEL 



We consider a fully-ionized quantum plasma composed of strongly correlated non-degenerate warm ion and mildly 
coupled degenerate electron fluids. Since, the ion thermal energy is much smaller than the electron Fermi-electron 
energy, we shall not consider quantum forces on the ion fluid. The viscoelastic hydrodynamic (VEHD) model for such 
an electron-ion plasma is governed by a set of equations Q 



(! + T vlk) [ m i n i^W ~ ZemE + VP; + Zn imi u ie Su\ 
= [r^VV + (£ + 2§l) V(V • uO] , Su = u t - u e 
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where r^, v e i{ v ie)i rj v , and Z are the viscoelastic relaxation time, electron- ion (ion-electron) collision frequency, 
shear viscosity, and the ion charge state, respectively. Also, the total electron pressure, Ptot = Pe + Pc + Pxo consists 
of pressures due to the electron degeneracy, Coulomb interaction, and electron-exchange interaction, respectively. The 
electron viscosity has been shown to become dominant only at much larger densities, such as in neutron-star crusts 
[30| . hence, we will ignore it in the following analysis. It is also known that in quantum plasmas the electron-electron 
collisions are rare due to the Pauli blocking-mechanism |5|. Therefore, it is convincing to only consider the ion-fluid 
viscosity effects in our model. 

We simplify the model by supposing that degenerate electrons are inertialess and the plasma is quasi-neutral, viz. 
n e ~ Zrii. Thus, we have p ~ mini, and obtain 



0, 



VP G ] = [ % V 2 u+(£+f)V(V-u)] 



(2) 



where, Pq — Ptot + Pi is the generalized plasma pressure. Moreover, it is convenient to express all involved 
parameters in terms of the Chandrasekhar's relativity parameter, R = PF e /m e c = (p/pc) 1 ^ 3 — (pe/ Pe) 1 ^ 3 , where 
p c ~ 2 x 10 6 gr/cm 3 and \x e = A/Z ~ 2 are the normalizing plasma mass-density and number of nucleons per electron 
with A and Z being the atomic weight and number, and Pp e is the electron relativistic Fermi momentum. 

Other quantities can be expressed in terms of the relativity parameter R and the plasma atomic-number Z. For 
instance, for the electron degeneracy pressure, we have 31 1 
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while, for the Coulomb interaction pressure in the spherical Wigner-Seitz cell approximation, we have [3S 
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Finally, for the electron exchange interaction [32| , one may write 
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where, j3 = R + yl + R 2 . On the other hand, the ion- viscosity is [33| 
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FIG. 1. Figure 1 shows the variations in strongly-coupled viscosity parameters of relativistically degenerate plasma ion- 
fluid, such as the viscoelastic relaxation-time, r, and the plasma ion fluid-viscosity, 77 with respect to different mass-density 
(relativistic-degeneracy parameter, Ro) and the atomic-number, Z. The curve thickness in plots indicate the increase in the 
varied parameter in each plot. 
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where, T = (Ze) 2 / '(roksT) is the coupling parameter with r = (3/47™ ) 1 / 3 being the Wigner-Seitz spherical-cell 
radius containing only one electron on average. Variations of the parameters I\ and I 2 are given in Ref. [33j for various 
coupling parameter values. According to DeWitt (3^ . the strongly coupled plasma is in liquid form for the range of 
the plasma ion coupling of Y < 178 to which we focus our attention in the following analysis. The plasma coupling 
factor, in terms of the relativity parameter for a typical white dwarf temperature of T ~ lO 6 ^ is approximately 
r = 22.752' 5 / 3 i?, valid for a wide ange of the plasma mass-density and atomic-number variations. The viscoelastic 
relaxation-time is directly related to the plasma viscosity via the following relation j33| 
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where, m* is the free electron effective-mass. The viscoelastic relaxation time is thus given as the following, with 
the fitting constants presented in Ref. [33| 
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Furthermore, ignoring the ion pressure compared to that of the electrons, we may write the plasma generalized 



pressure in terms of a generalized effective potential 35 1 
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with the effective plasma potential, defined as 
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where, /3 = 3 1/3 (2a/5)(2Z/7r) 2 / 3 . The dimensionless set of VEHD equations, Eqs. ©, can be obtained by using the 
following scalings 

c - i 

V -> — V, t-> , p^ppo, u^uc s , R->Rqp 1/3 , (11) 

where, the normalizing factors, po, Lu P i — (e/m^^/Anp^ ~ 10 18 Z\/i?o (with i?o = (po/Pc) 1 ^ 3 and p c ~ 2 x 10 6 <7r/cm 3 ), 
and c s = cy m e /mj, denote the equilibrium plasma mass-density, characteristic ion plasma frequency and the plasma 
sound speed, respectively. Thus, the normalized equations, dropping the bar notation for simplicity, read 



d t p + pV ■ u = 0, , , 

(1 + rdt) [p (dtu + W G )] = r, [V 2 u + |V(V ■ u)] . [iZ) 

where, r = a; pi T„ ~ 1.92y/B$/(I 2 y/l + R%f? /3 ) and ^ = (w^/c^ ~ 2.07i?o /2 /[^/5 1/3 (l + i?gp 2/3 )]. It should 
be noted that, at very high densities, such as in the neutron star crust region, the electron-electron collisions can 
be very important dominating the electron viscosity over that of the ions, the consideration of which is beyond the 
scope of the present investigation. Equations (|12[) can be considered as the most general model for an unmagnetized 
viscoelastic fully interacting quantum plasma, valid for a wide- range of the plasma mass density, atomic-number, 
and the ion-fluid coupling factor. The quantum plasma parameters range of interest covering the planet interiors and 
white dwarfs are approximately, 0<i?o<l,l<^<12 and 1 < V < 178. The quantum fluid model, following 
the closed equation-set relevant to different plasma geometries for the one-dimensional (planar (y = 0), cylindrical 
(u = 1), and spherical {v = 2)) coordinate system, can be written as 



dp(r,t) 

at 



[l+T(p,Z)£] [p(r,t) + 

i\r-> i r v Qr \ or J 3 or \r u or J 



(13) 



In the following sections, we will evaluate the characteristics of the viscoelastic ES nonlinear excitations in terms of 
fractional Fermi-Dirac plasma parameters. Note that, in our model, we have ignored the bulk viscosity, C> compared 
to that of shear, 77, due to negligible contributions from compressibility of high pressure plasma fluids. 

Figure 1 displays the variations of (a) the normalized viscoelastic relaxation time, r, and (b) the normalized shear 
ion fluid- viscosity, 77 with respect to the plasma degeneracy parameter, Rq, for different plasma atomic-number, Z. 
It is observed from Fig. 1(a) that for all atomic- number values the viscoelastic relaxation time increases with the 
increase of the plasma number density, related to the relativistic degeneracy parameter Rq, and this increase is smaller 
for larger values of plasma atomic- number, Z. Figure 1(b), on the other hand, reveals that the plasma ion- viscosity 
first increases with the increase of the plasma number density up to a maximum value and then decreases again. 
This behavior is observed for all plasma atomic-number value being more pronounced for plasmas with lower atomic 
numbers. It is also remarked that the relative value of the ion fluid- viscosity parameter for a given plasma mass 
density is always lower for heavier element composed plasmas. 



LINEAR AND NONLINEAR DENSITY PERTURBATIONS 



Important information can be obtained considering the plasma response to linear perturbations. This is accom- 
plished by the Fourier analysis of Eq. (|13j) . and by using the harmonic operators d r — ik and dt — —iui. Accordingly, 
we have the following first order relations 



kui — topi = 0, 

3(t + TLj)(kpiT - uiuj) + irjk 2 ui = 0, 



(14) 
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FIG. 2. Figure 2 shows the linear wave dispersion of ion waves in the relativistically degenerate strongly-coupled viscous plasma 
and the variation of wave speeds in terms of the ion mass-density (relativistic-degeneracy parameter Ro) and the atomic- number 
Z. The curve thickness in plots indicates the increase in the varied parameter in each plot. 
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From Eqs. (|14|l we readily obtain the generalized dispersion relation 



3ik 2 T + k 2 (4?7 + 3Tt)u> - iiuj 2 - 3rw 3 = 0. (16) 

Figure 2 depicts the linear dispersion of viscoelastic ion waves in a relativistically degenerate plasma and the 
variation of the phase-speed for different plasma mass density and atomic-numbers. Figure 2(a) reveals that the 
linear wave dispersion is fundamentally governed by the plasma mass density. Below a critical plasma density, a 
dispersion gap opens in the wave-number range, indicated for the thin curve in the plot. As the plasma density 
increases above the critical value, for a fixed plasma composition, the dispersion gap disappears, but there still exists 
a sharp discontinuity in the linear wave speed. Also, as the plasma mass-density increases further this discontinuity 
disappears. Such anomalous feature may be attributed to the interplay of strong plasma interactions and ion fluid 
couplings expected for a super-dense degenerate plasma. Figure 2(b) shows the effect of the plasma composition on 
the linear wave dispersion. It is evident that in the small and long-wavelength limits, the linear wave speed is not 
affected by the change to plasma composition, while this is converse in the medium-wavelength region due to the 
presence of the mentioned discontinuity in the plasma wave speed. 

It is evident that the generalized set of VEHD equations, Eqs. (TT5T) . do not admit analytical arbitrary amplitude 
solutions. Therefore, at this stage, we aim at numerical investigation of arbitrarily large amplitude nonlinear waves 
which can have fundamental importance in superdense astrophysical environments. In order to evaluate the nonlinear 
blast- waves in the laboratory frame, we solve the time-dependent VEHD equations compliant with appropriate periodic 
boundary conditions for different independent fractional plasma parameter values, namely, the relativistic degeneracy 
parameter i?o (representing the plasma mass density, p), and the plasma composition parameter, Z. For the the 
planar geometry (viz. v = 0) a Gaussian pulses of the form p(r,0) = 1 + exp(— r 2 /100) with an initial speed profile 
of u{r, 0) = 0.7exp(— r 2 /100) is placed at the origin in initial time. The time evolution of the plasma mass-density 
perturbation is then rendered to future times based on the governing differential equation set, Eq. (|13[) . for given 



7 



R = 0.9, Z = 2, r = ?7 = 








FIG. 3. Figure 3 shows the evolution of a Gaussian density pulse in a non-viscous (r? = r = 0) collisional quantum plasma for 
a given plasma number density (no ~ 4.3 x 10 29 cm -3 ) and atomic-number (Z = 2). 



R = 1, Z = 2, T = 10 6 (K) /?o = 1, Z = 2, T = 10 6 (K) 




FIG. 4. Figure 4 shows the evolution of a Gaussian density pulse and the corresponding ion-speed profiles in a strongly-coupled 
collisional quantum plasma (with a given ion mass density, atomic-number and the plasma temperature) into a blast-like 
nonlinear ion- wave structure. 

values of the plasma fractional parameters. It is further assumed that the spatial variations in the plasma mass-density 
vanishes at the plasma boundaries in all future times. In our model, we also incorporate the effects of time-dependence 
of the viscoelastic parameters due to the plasma mass-density variation, for completeness. The results of simulation 
are depicted in Figs. 3 and 4 for non-viscous (77 = r = 0) and viscous (t; ^ t ^ 0) quantum plasmas, respectively, 
for given initial plasma mass-density, temperature and composition relevant to typical white dwarf stars and cores of 
giant planets. Figure 3 indicates the time evolution of ion density wave into the train of density humps with the time 
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FIG. 5. Figure 5 shows the evolution of a presummed Gaussian density profile in a strongly-coupled collisional quantum plasma 
for different simulation times and for different plasma fractional parameters. The figure consists of different plots for variations 
of the plasma fractional parameters, namely, the ion mass-density (denoted by the relativistic-degeneracy parameter Ro) and 
the atomic-number Z. Each row in the figure indicatess the effect of the variation in atomic number on the shock height and its 
steepness and each column depicts the effect of the change in the ion mass density on the shock wave characteristic parameters. 
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R = 1, Z = 2, T = 10 5 (K), v = R = 1, Z = 2, T = 10 5 (K), v = 1 




FIG. 6. Figure 6 shows the existence of large-amplitude nonlinear electrostatic shock excitations in different geometries, namely, 
planar (y — 0), cylindrical (i> = 1), and spherical (i/ = 2). Plot (d) reveals that the shock height increases as one goes from 
spherical to cylindrical and from cylindrical to planar geometry. 



elapse. This is an example of ion wave breaking/steepening, arising from harmonic generation nonlinearity, where, 
the plasma dispersive effects are somewhat spreading the ion pulse, but not strong enough to allow the formation of 
a real solitary pulse. These ion density pulses arise from balance between harmonic nonlinearity and the quantum 
statistical pressure and interaction effects in the limit of vanishing plasma viscosity and the corresponding relaxation 
time. For the case of strongly coupled quantum plasma the value of the plasma coupling factor, T, is kept below 178 
in our simulations to ensure that the quantum plasma is in the liquid phase 34]. It is observed from Fig. 4 that the 
time evolution of the Gaussian density and velocity profiles develop into shock-like structures with pronounced shock 
wave front steepened by further elapse of the simulation time. 

It is confirmed by Fig. 5 that for all the considered set of plasma parameters the Gaussian perturbation profile 
develops into shock-like structures after enough elapsed simulation time, say (t ~ 100/cj p ,). Comparing different plots 
in Fig. 5, we observe that an increase in both the plasma mass density and the atomic-number leads to an increase in 
the shock height and steepening of the shock front. Remember that the maximum value of the relativistic degeneracy 
parameter, Rq — 1 used in this simulation, corresponds to high electron/ion number density, n ~ 5 x 10 29 cm -3 , 
expected for white dwarf stars. The present simulation revealed that despite the significant effect of the temperature 
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on the plasma coupling factor, T, it has little effect on the viscosity parameters and the present study might be 
extended to the (Jupiter-like) large planet cores and environments with high number densities and relatively low 
temperatures. 

Figure 6 depicts the evolution of ion-momentum distribution from an initial Gaussian profile in different plasma 
geometries. It is seen that in planar [y = 0), cylindrical (u = 1) and spherical [y — 2) geometries, there exist 
large amplitude shock waves. In producing Fig. 6(a-c), we have assumed an initial pulse shape as p(r, 0) = 1 + 
exp[— (r — 100) 2 /100] and u(r, 0) = 0.7exp[— (r — 100) 2 /100] in order to avoid singularities present at the origin for 
non-planar geometries. From Fig. 6(d) we observe that the amplitude of nonlinear ion shock structures are relatively 
higher in the planar geometry with respect to other geometries, for a given set of fixed plasma parameters. Different 
geometries than planar, like spherical and cylindrical, may be important in the study of physical mechanisms in 
stellar configurations and active galactic nuclei, for instance. The present findings reveal the fundamental effects of 
the plasma mass density and its composition on the formation and characteristics of blast-like ion-waves in dense 
quantum plasmas. Furthermore, our study reveals that the characteristics of density perturbations and the formation 
of ES ion shocks in the cores of white dwarf stars can be quite different from the existing shock-waves at the low 
density crust region. 

SUMMARY AND CONCLUDING REMARKS 

In this paper, we have presented a nonlinear theory for arbitrary large amplitude electrostatic planar and non- 
planar shocks in collisional quantum plasmas composed of mildly coupled degenerate electron fluid and strongly 
correlated non-degenerate ions for a wide range of the plasma mass density and atomic number. We have included 
the electrostatic and relevant pressure gradient, quantum forces and interactions on degenerate electron fluids, and 
have used the general VEHD equations for non-degenerate ions to carry out laboratory-frame simulations in different 
geometries in order to deduce the complete picture of very large-amplitude electrostatic shock structures. The latter 
may exist in white-dwarfs, the active galactic nuclei or many other high density astrophysical environments such 
as neutron stars, pulsars and magnetars. It was found that that the increase/decrease in the ion charge state or 
atomic- number/mass-density would lead to stronger and steeper shock fronts in the relativistically degenerate viscous 
plasmas. It was also shown that the linear dispersion in such quantum plasmas exhibits anomalous features attributed 
to the strong fermionic interactions and ion fluid couplings. In conclusion, we stress that the present investigation of 
the formation of arbitrary large amplitude ES ion shocks will play a very important role in acceleration of electrons 
and ions in relativistically degenerate quantum plasmas, which are ubiquitous in high-energy density compressed 
plasmas produced by intense short wavelength laser beams, as well as in planetary systems (e.g. the giant Jupiter) 
and in super dense compact astrophysical objects (e.g. white dwarf stars and magnetars). 
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